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(1. 华北 电力 大 学 核 科 学 与 工程 学 院 北京 102206) 
摘要 : 为 提高 空间 堆 中 子 学 模拟 的 可 靠 性 ， 基 于 国际 临界 安全 手册 中 关于 空间 堆 


的 临界 基准 实验 (HMF005 基准 实验 ) ， 采 用 蒙特 卡 罗 程 序 对 HMF005 基准 实验 
中 堆 芯 排 布 进 行 建 模 ， 结 合 不 同 评价 核 数 据 库 ， 计 算 6 种 堆 芯 排 布下 的 kr 计算 
值 与 实验 值 的 偏差 ， 并 对 计算 值 与 实验 值 的 比值 的 平均 值 、 标 准 偏差 、 交 和 均 方 
根 误差 等 进行 深入 分 析 。 研 究 发 现 ，ENDF/B-VIL1、CENDL-3.2、ENDF/B-VIIL0 
和 JEFF-3.3 四 种 评价 核 数 据 库 计算 扒 芯 kg 结果 与 实验 值 偏差 处 于 三 倍 不 确定 度 
之 内 ， 其 中 ENDF/B-VL1 计算 结果 最 接近 实验 值 ，CENDL-3.2 与 实验 偏差 稍 大 。 
CENDL-3.2 库 计算 结果 与 实验 值 偏 差 稍 大 主要 是 由 于 评价 库 中 的 钼 和 皱 核 素 导 
致 ， 钼 的 同位 素 中 "Mo 的 影响 最 大 。 同 时 还 进行 敏感 性 分 析 ，%Mo 和 ?Be 的 总 
上 反应 敏感 性 系数 主要 跟 弹 性 散射 敏感 性 系数 相关 , CENDL-3.2 和 ENDF/B-VIL.1 中 
Be 的 辐射 俘获 敏感 性 差异 较 大 ， 两 者 偏差 可 达 34.995. 
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中 图 分 类 号 : TL375.5 ”文献 标志 码 : A doi: 


EH 
Till 


随 着 世界 航天 探索 技术 的 发 展 , 目前 , 基于 太阳 能 、 化 学 能 的 空间 电源 技术 能 力 已 经 接 
近 极 限 水 平 , 随 着 空间 活动 不 断 深 入 ,常规 能 源 将 难以 满足 深 空 探测 和 特殊 任务 的 发 展 要 求 ; 
而 核能 因 其 具有 功率 覆盖 范围 广 、 能 量 密度 高 、 受 环境 影响 小 、 寿 命 长 等 特点 ,在 空间 探索 
中 具有 广泛 的 应 用 前 景 。 目 前 , 美国 、 俄 罗斯 和 欧盟 等 国家 都 把 空间 核 动力 技术 视 为 国家 前 
沿 核心 技术 ， 并 积极 围绕 空间 应 用 和 空间 资源 开发 进行 战略 部 署 03。 由 于 利用 数字 反应 堆 
技术 可 以 缩短 反应 堆 研 发 周期 和 降低 研发 成 本 ,目前 各 国都 在 利用 数字 反应 堆 技 术 进 行 空间 
堆 设计 , 评价 核 数据 库 中 数据 的 精度 则 直接 影响 空间 反应 堆 设 计 精 度 , 因此， 开展 评价 核 数 
据 库 应 用 于 空间 核反应 堆 设计 精度 的 评估 对 新 型 空间 堆 的 设计 具有 重要 意义 ,虽然 国 内 外 对 
空间 堆 有 大 量 研究 ， 比 如 前 苏联 在 60 年 代 陆 续 研 制 出 TOPAZ-I 和 TOPAZ-I 等 型 号 的 空间 
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IZMED, 美国 在 2018 年 完成 千瓦 级 空间 堆 Kilopower 反应 堆 的 测试 内， 中 国 原 子 能 院 提 
出 SPACE-R ZEHE, 清华 大 学 REAL 团队 开发 的 ACMIR 型 空间 堆 回 , 北京 航天 器 环境 工 
程 研究 所 可 靠 性 与 环境 科学 技术 实验 室 研发 的 ACMIR 型 空间 堆 门 ， 南 京 航空 航天 大 学 的 李 
婷 等 人 提出 的 熔 盐 冷却 空间 堆 概 念 四 等 。 除 了 之 前 的 实验 研究 ， 目 前 很 多 的 概念 设计 大 多 直 
接 使 用 评价 核 数据 库 进 行 空间 堆 反 应 堆 物 理 计算 ,很 少 对 评价 核 数据 库 在 空间 堆 设 计 中 的 可 
靠 性 进行 评价 。 因此 开展 空间 堆 核 数据 宏观 检验 研究 对 空间 反应 堆 的 设计 非常 重要 。 临界 基 
准 实验 通常 被 用 于 评价 物理 计算 模型 、 设 计 方法 以 及 所 用 的 核 数 据 库 的 准确 性 ; 基于 实验 精 
确 测量 反应 堆 的 re 可 实现 对 堆 臣 关键 材料 核 数 据 的 宏观 基准 检验 外。 为 了 对 空间 堆 反 应 堆 
物理 计算 使 用 的 评价 核 数 据 库 进行 临界 验证 ， 本 文通 过 选取 国际 临界 安全 手册 ( ICSBEP ) 
中 HMF005 HEERO, RH MCNP 程序 以 及 不 同 评价 核 数据 库 ， 计 算出 与 实验 对 应 的 堆 世 
kr， 开展 空间 堆 相 关 的 核 数据 入 堆 宏 观 检 验 ， 通 过 分 析 不 同 核 数据 库 ke 计算 值 与 实验 值 差 
值 的 大 小 ， 检 验 目前 评价 核 数据 库 的 可 靠 性 。 利 用 了 MCNP 程序 计算 钼 和 钙 的 敏感 性 系数 
和 不 同 堆 芯 模型 下 的 中 子 通 量 密度 。MCNP 使 用 了 伴随 加 权 的 线性 扰动 理论 计算 敏感 性 系 
1 临 青 基准 实验 

HMF005 系列 实验 为 1987 年 俄罗斯 奥 布 守 斯 克 物 理 和 动力 工程 研究 所 进行 的 临界 实验 。 
实验 按照 堆 芯 材料 排 布 不 同 , 可 以 分 为 六 种 堆 必 。 六 种 反应 堆 堆 世 都 采用 高 浓缩 铀 金属 作为 
燃料 ， 由 于 铀 合金 具有 高 反射 率 、 耐 高 温 、 耐 腐蚀 等 优点 ， 部 分 反应 堆 设 计 采 用 金属 钼 反射 
JZ RTE, 部 分 反应 堆 设 计 也 采用 争 作 为 反射 层 。 以 上 反应 堆 设 计 也 为 检验 钼 和 争 的 评价 
中 子 截面 数据 提供 很 好 的 实验 数据 HMF005 系列 基准 题 包括 六 个 实验 , 本 文 使 用 蒙特 卡 罗 
程序 对 六 个 实验 模型 进行 精细 化 建 模 。 表 1 中 给 出 组 成 实验 中 用 到 各 类 棱 束 的 基本 组 成 材 
料 以 及 每 种 材料 的 原子 密度 , 实验 中 共 涉 及 到 5 种 材料 , 分 别 是 铀 块 、 反 射 层 ( 底部 和 顶部 )、 
THEE. IER. MEE, CHE RHEL IM U3 HR. U6 PE. US FE, UII HE, THEE. BOR ME 
管 等 不 同类 型 的 燃料 或 者 其 它 棒 。 SEES Ua B pip ZH Ah. 05] SERIES TEE HER 
六 个 堆 芯 实验 模型 ， 具 体 堆 芯 排 布 和 堆 芯 组 成 如 下 图 1- 图 6 所 示 ， 尺 寸 单位 为 厘米 。 

表 1 各 区 域 材料 原子 密度 (1/barn'cm) 


[ES 铀 块 顶部 和 低 部 反射 层 钼 棒 pi 空 管 
a 3.3109x104 
SU 3.2968x102 
as 1.4590x 10^ 


228Ly 3.1467x10? 


5.3547x10? 5.3547x10? 


3.0875x10? 


Molybdenum 


1.0962x 107 


Beryllium 


5.3293x10? 3.2273x10? 


3.2273x10? 5.3293x10? 


3.2912x10? 


Tron 


8.4596x10^ 


1.3969x10? 


8.4596x10^ 8.4596x10^ 1.3969x10? 


Chromium 


7.5543x 104 4.5748x104 


4.5748x 104 4.5748x 104 7.5543x104 


Nickel 


9.2673x 10% 5.6121x105 


5.6121x105 5.6121x105 9.2673x105 


Manganese 


7.5532x105 4.5741x105 


4.5741x105 4.5741x105 7.5532x10° 


Silicon 


4.2090x 10° 2.5489x 10° 


2.5489 105 2.5489x 10° 4.2090x 10° 


Titanium 


3.7090x105 2.2461x105 


1.0624x10? 2.2461x105 3.7090* 10° 


Carbon 


3.0863x10? 


Aluminum 
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图 1 HMF005_01 堆 芯 排 布 ( 左 ) 和 堆 芯 组 成 ( 右 ) 
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HMFO005. 02 堆 世 排 布 ( 左 ) 和 堆 芯 组 成 ( 右 ) 
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图 3 HMF005_03 HEHE (Æ) 和 堆 芯 组 成 〈 右 ) 
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4 HMFO005. 04 推 芯 排 布 〈 左 ) 和 推 芯 组 成 ( 右 ) 
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图 5 HMFO005. 05 HERES (Ze) 和 堆 芯 组 成 ( 右 ) 
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图 6 HMFO005. 06 堆 芯 排 布 〈 左 ) 和 堆 芯 组 成 ( 右 ) 


本 文 首 先 基于 评价 核 数据 库 ENDF/B-VI.1!! , CENDL-32U2 | ENDE/B- VII.0U3! 和 
JEFF-3.304 利用 评价 核 数 据 库 处 理 软件 NJOY2016 制作 了 MCNP 程序 使 用 的 连 


据 库 ， 再 使 用 MCNP 程序 对 有 效 增殖 系数 Cs 以 及 相关 核 素 的 敏感 性 


点 截面 数 


4 


续 


系数 进行 计算 。 


准 偏差 .计算 公式 如 下 : 


AN 


评价 参数 选用 计算 值 与 实验 值 的 比值 的 平均 值 以 及 其 


<C/E>=— - YE (1) 
i=l 
ac = Lt CIE <C/E>)’ (2) 
n-— i-l 
Hop, C WE , 分 别 代 表 对 应 实验 i 的 ks 计算 值 和 实验 值 。 同 时 ，X2 和 均 方 根 误差 


(RMSE ) 也 用 于 评估 计算 结果 和 实验 


RB AE, 表示 计算 值 与 实验 值 的 偏 
模拟 计算 时 设 
代 ， 非 活跃 代为 100 代 。 
3 计算 结果 分 析 
3. 1 有 效 增殖 系数 


0.990 F 


0. 985 


SS WUE A 300K. FEI ARITA, UBS10000 个 粒子 ， 


吉 果 之 间 的 一 致 性 ， 具 体 表示 如 下 : 


ICO Ea 
= i i 3 
aD AES ) (3) 
RMSE = Ey Qn? (4) 
n 
AE, = C, — E, (5) 


差 。 基 准 实验 是 在 300K 温度 下 进行 , 故 本 文 使 用 MCNP 
总 代数 为 1100 


ICSBEP-2014 ”实验 值 
ENDF/B- vi 0 n 算 信 
CENDL-3 KACI 
ENDF/B-VILI inet 
JEFF-3.3 算 值 


e4pen 


图 7 24 


ak 2 HMF005 基准 题 计算 结果 的 x 


F 不 同 评价 核 数据 库 计 算 值 与 实验 值 对 比 


2 , RMSE, <C/E>-1 以 及 A (C/E) 


ENDF/B-Ml.1 CENDL-3.2 ENDF/B- VII.O JEFF-3.3 
x 2 0.686553498 3.945609568 0.972362397 1.850668724 
RMSE 298.290686 715.0881065 354.9903755 489.7414284 
C/E -1 -0.001903728 -0.006670100 -0.002709571 -0.004204959 
A(C/E) 0.002517317 0.005940129 0.002666283 0.002753523 


图 7 给 出 基于 不 同 评价 核 数 据 库 采用 MCNP 程序 对 六 个 实验 计算 的 kr 结果 ， 其 中 基 
准 题 实验 结果 的 不 确定 度 为 360 pem， 六 个 基准 题 蒙 卡 计算 结果 的 统计 误差 最 大 不 超过 21 
pem. 由 图 7 AIL, 部 分 评价 核 数 据 库 的 计算 结果 在 1 倍 置 信 区 间 之 内 ,所 有 计算 结果 都 处 
于 3 倍 置信 区 间 之 内 。 根 据 第 二 节 中 的 分 析 方 法 ， 将 采用 不 同 评价 核 数据 库 计 算得 到 的 6 
个 堆 芯 的 kr 结果 与 基准 题 实验 值 进行 处 理 ， 结 果 如 表 2 所 示 。 由 表 2 可 见 ， 采 用 
ENDF/B- VI.1 库 和 ENDF/B-VILO 库 的 计算 结果 与 实验 值 结果 符合 较 好 ,但 基于 CENDL-3.2 
评价 核 数据 库 的 结果 明显 与 其 它 三 个 库 结果 均 有 较 大 差别 ， 均 方 根 误差 (RMSE ) 大 约 是 
ENDF/B-VIILO 计算 结果 的 2 倍 。 由 于 基准 实验 中 大 量 采用 钼 和 争 核 素 做 反射 材料 ， 上 述 差 
异 可 能 是 由 于 反射 层 材料 评价 核 数据 的 不 同 导致。 为 研究 CENDL-3.2 库 结果 偏差 较 大 的 原 
, 将 基准 实验 中 6 种 堆 芯 排 布 的 反射 材料 全 部 去 掉 , 再 使 用 CENDL-3.2 PERI ENDF/B- VI.1 
库 进行 Cr 计算 ， 具 体 实 验 结果 如 表 3 所 示 。 由 表 3 可 见 ， 去 掉 反 射 棒 以 后 ，CENDL-3.2 
库 与 ENDF/B-W.1 库 的 计算 吻合 非常 好 ， 由 此 可 以 推断 ，CENDL-3.2 库 与 其 它 评价 核 数据 
库 的 偏差 主要 自 于 反射 层 材 料 钼 棒 和 钙 棒 评价 核 数据 的 不 同 。 


表 3 HMF005 基准 题 无 反射 层 CENDL-3. 2 和 ENDF/B- VII. 1 k,,, 计算 结果 


01 模型 ”02 模型 03 模型 O4 模型 05 模型 ”06 模型 


CENDL-3.2 0.73984 0.69527 0.66558 0.64135 0.69282 0.71136 
ENDF/B-VI.1 0.73966 ^ 0.69557 0.66535 0.64153 0.69564 0.71098 
Kr 偏差 /pcem — -18 30 23 18 -18 -38 


在 HMF005 04 HAA 1 根 钼 反射 棒 ， 可 视 HMF005 04 HERET > ESI 
HMFO005 01 (全 钼 棒 反 射 ) 和 HMF005 04 ( 2 BRED) 的 特殊 性 ， 因 此 对 这 两 个 基准 实 
验 进行 深入 研究 。 以 CEDNL-3.2 计算 结果 为 基础 ， 用 ENDF/B-VI.1 库 中 的 钥 和 皱 蔡 换 掉 
CENDL-3.2 PE PAGER, SEAT HMF005_01 基准 题 和 HMF005_04 基准 题 计算 ， 
计算 结果 与 实验 值 的 偏差 对 比 结果 如 表 4 所 示 。 由 表 4 可 见 ， 将 HMF005_01 中 的 钥 核 素 
和 HMF005_04 rp AU BlEgz 3521 ERHÁ. y ENDF/B- VI.1 库 来 源 后 ， 使 用 CENDL-3.2 库 计算 的 
kf 与 实验 值 的 偏差 明显 偏 小 ， kr 分 别 改 进 643 pem 和 402 pcem。 通过 控制 单一 变量 的 原则 ， 
对 HMF005_01 中 钥 核 素 的 天 然 同 位 素 分 别 进 行 核 数据 库 的 替换 来 研究 钼 各 同位 素 对 ker 
来 的 影响 , 计算 结果 图 8 所 示 。 由 图 8 可 见 , 对 于 HMF005_01 基准 题 , 钼 的 同位 素 中 , ?Mo 
的 影响 最 大 ， 为 301 pom; Mo 的 影响 最 小 ， 其 它 同 位 素 大 约 有 100 pem 的 影响 。 

表 4 计算 值 与 实验 值 ki 的 偏差 (C-E) 对 比分 析 (单位 : 1. 0X10”) 


es CENDL-3.2( 反 射 层 中 的 钼 核 数 据 (HMF005_01) 或 者 锌 数据 
评级 核 数据 库 ENDF/B-VI.1 CENDL-3.2 CHMF005_ 04) 采用 ENDF/B-VII.1 中 的 相应 数据 蔡 换 ) 


HMF005501 457 1038 395 


HME005_04 473 919 517 
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图 8 HMF005_01Mo 各 同位 素 引 起 的 CENDL-3.2 与 ENDF/B-VI.14,, 25 fH 


3. 2 敏感 性 系数 分 析 

针对 HMF005 01 和 HMF005 04 数据 库 替 换 后 带 来 kg 的 差异 ， 对 ”Mo Fil °Be 在 
ENDF/B- VI.1 库 和 CENDL-3.2 库 的 总 反应 、 弹 性 散射 、 总 吸收 和 辐射 俘获 反应 敏感 性 系 
数 进行 研究 。HMF005 01 中 ?Mo 四 种 反应 下 的 总 敏感 性 系数 如 表 5 所 示 ， 对 比 发 现 
CENDL-3.2 和 ENDF/B- VI.1 中 ?Mo 的 总 反应 敏感 性 系数 和 总 吸收 截面 的 总 敏感 性 系数 
差异 都 在 1% 以 内 ,弹性 散射 系数 和 辐射 俘获 的 总 敏感 性 系数 偏差 达到 7%。 从 图 9- 图 12 
可 知 ，%Mo 总 反应 敏感 性 系数 主要 由 弹性 散射 敏感 性 系数 决定 ，”2”Mo 核 素 辐射 俘获 和 
吸收 截面 的 敏感 性 系数 较 小 , 在 中 高 能 区 ( 中 子 能 量 大 于 104 ev ) , Mo 在 ENDF/B- VI.1 
和 CENDL-3.2 敏感 性 系数 有 显著 差异 ， 这 可 能 导致 kr 差 值 较 大 的 原因 。 


表 5 CENDL-3. 2 和 ENDF/B- VII. 1"Mo 四 种 反应 道 总 的 敏感 性 系数 


评价 核 数 据 库 ”总 反应 敏感 性 ”弹性 散射 敏感 性 ”吸收 敏感 性 ”辐射 俘获 敏感 性 


CENDL-3.2 3.418 X 10? 3.203 X 10? -1.581 X103 . -L506x10? 
ENDF/B-VLI ^ 3421X10? 3.237 X 107 -1.600X10?  -1.612X10? 


相对 误差 0.109% 1.04% 6.48% 6.57% 
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图 11 Mo 总 吸收 敏感 性 系数 12 Mo 辐射 俘获 反应 敏感 性 系数 


HMF005_04 中 ”Be 在 四 种 反应 下 的 敏感 性 系数 如 表 6 所 示 ， 对 比 发 现 CENDL-3.2 和 
ENDF/B- VI.1 中 ?Be 的 总 敏感 性 系数 、 弹 性 散射 系数 和 总 吸收 系数 差异 都 在 2% 以 内 ,但 辐 
射 俘获 敏感 性 系数 差别 较 大 ; 图 13- 图 16 分 别 给 出 了 总 反应 截面 、 弹 性 散射 截面 、 吸 收藏 
面 和 辐射 俘获 截面 的 敏感 性 系数 分 布 。 由 图 13- 图 16 ATL, Be 的 总 反应 敏感 性 系数 主要 
跟 弹 性 散射 敏感 性 系数 相关 。CENDL-3.2 和 ENDE/B- VI.I 中 "Be 的 辐射 俘获 敏感 性 差异 较 
大 ， 两 者 偏差 可 达 34.9%。 虽 然 两 个 评价 核 数据 库 "Be 的 辐射 俘获 敏感 性 差异 较 大 , 但 由 于 
辐射 俘获 反应 道 的 敏感 性 系数 较 小 , 占 吸收 截面 的 敏感 性 系数 较 小 , 因此 ,吸收 截面 的 敏感 
性 系数 差别 不 大 。 


表 ó CENDL-3. 2 和 ENDF/B- VII. 1°Be 四 种 反应 的 总 的 敏感 性 系数 


评价 核 数据 库 ”总 反应 敏感 性 弹性 散射 敏感 性 吸收 敏感 性 辐射 俘获 敏感 性 
CENDL-3.2 0.2878 0.267 -4.927 X 10° -2.644 X 104 
ENDF/B-VIl.1 0.293 0.272 4.840 X 103 -1.960 X 104 


相对 误差 1.80% 1.61% 1.81% 34.9% 
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图 15 "Be 总 吸收 敏感 性 系数 16 "Be 辐射 俘获 反应 敏感 性 系数 

3.3 归 一 化 中 子 通 量 密度 

基于 不 同 评价 核 数 据 库 对 六 个 基准 实验 进行 全 堆 世 的 中 子 通 量 密度 计算 ,计算 结果 如 下 
图 17- 图 22 所 示 。 中 子 通 量 密度 采用 30 群 能 群 结构 ， 具 体 的 能 群 结构 可 以 参考 [10]。 由 图 
17- 图 22 可 见 , 基于 不 同 评价 核 数 据 库 ENDF/B-VI.1 , CENDL-3.2 , ENDF/B-VII.O 和 JEFF-3.3 
计算 的 六 种 堆 芯 下 的 归 一 化 中 子 通 量 密度 趋势 一 致 ， 差 别 较 小 。 由 HMF005_02 基准 实验 到 
HMF005_06 基准 实验 可 以 明显 看 出 反射 层 Be 对 中 子 的 慢 化 作用 ,而 模型 1 相对 于 其 它 基准 
模型 ， 中 子 能 谱 更 硬 。 
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图 17 HMF005.01 归 一 化 中 子 通 量 密度 18 HMF005_02 归 一 化 中 子 通 量 密度 
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图 21 HMF005 05 归 一 化 中 子 通 量 密度 ”图 22 HMF005_06 归 一 化 中 子 通 量 密度 


4 结论 
基于 国际 临界 安全 基准 手册 基准 题 (HMF005 实验 ) , 开展 了 空间 反应 堆 评 价 核 数据 的 

宏观 检验 研究 。 利 用 MCNP 程序 和 四 种 评价 核 数 据 库 ( ENDF/B-W.1 、CENDL-3.2、 
ENDF/B-WL0 和 JEFF-3.3 ) 计算 了 六 种 堆 芯 排 布 空间 反应 堆 基 准 实验 的 kr。 主 要 结论 为 : 

(1 ) 基准 题 的 实验 值 与 蒙 卡 程序 结果 显示 ， 蒙 卡 计算 结果 与 实验 结果 偏差 均 处 于 三 倍 标准 
偏差 之 内 。 其 中 ENDF/B-VIL.1 库 与 实验 值 吻 合 最 好 , CENDL-3.2 库 与 实验 值 相差 偏 大 ; (2) 
CENDL-3.2 库 计 算 结 果 与 实验 值 偏 差 较 大 主要 来 自 于 反射 层 材料 钥 和 皱 核 数 据 。 如 果 采 用 
ENDF/B-VI.1 中 钼 和 外 的 核 数 据 蔡 换 基于 CENDL-3.2 评价 库 的 连续 点 截面 数据 ， 
HMF005 01 ( 全 钥 棒 反射 ) 和 HMF005_04 (全 皱 棒 反 射 ) 的 kr 分 别 可 以 改进 643 pem 和 
402 pcem。 对 钼 的 天 然 同 位 素 单独 研究 后 ， 发 现在 CENDL-3.2 库 和 ENDF/B-W.1 库 中 核 素 
Mo 的 改变 引起 的 kg 变化 最 大 ; (3) 通过 对 ENDF/B- VI.1 和 CENDL-3.2 中 ”Mo 和 ?Be 
的 敏感 性 系数 进行 研究 ， 发 现 ”Mo 和 ?Be 的 总 反应 敏感 性 系数 主要 跟 弹性 散射 敏感 性 系数 
相关 ,CENDL-3.2 和 ENDF/B- VII.1 中 "Be 的 辐射 俘获 敏感 性 差异 较 大 ,两 者 偏差 可 达 34.9%。 
从 六 个 基准 实验 的 中 子 能 谱 可 见 ，?Be 对 中 子 的 慢 化 效果 明显 。 


参考 文献 


[1] HU Jian, SU Jiqiang, LI Yanrui, et al. Industry and Technology Forum,2022,21(6):33-35(in Chinese). 
DOLI:10.3969/j.issn.1673-5641.2022.06.013. 

( 骨 健 , 宿 吉 强 , 李 言 瑞 ,等 。 产业 与 科技 论坛 ,2022,21(6):33-35. DOI10.3969/j.issn.1673-5641.2022.06.013.) 

[2] ZHOU Biao, JI Yu, SUN Jun, et al. Atomic Energy Science and Technology,2020,54(10):1912-1923(in Chinese). 
DOI:10.7538/yzk.2019.youxian.0736 


( 周 彪 , 吉 宇 , 孙 俊 ,等 . 原子 能 科学 技术 ,2020,54(10):1912-1923. DOI:10.7538/yzk.2019.youxian.0736.) 


[3] Susan S. Voss. Nuclear Technology,2022,206:8,1097-1108. 
[4] ZHANG Yan,WU Haosong, WANG Shu. National Defense Science and Technology Industry,2018(7):65-66(in Chinese). 


( 张 焰 , 伍 浩 松 , 王 树 . 国 防 科技 工业 ,2018(7):65-66.) 
[5] LIU Lili, SUN Zheng, FU Ziming. Atomic Energy Science and Technology,2015(9):1624-1628(in Chinese). 
DOI:10.7538/yzk.2015.49.09.1624. 


( 刘 黎 丽 , 孙 征 , 付 子 明 . 原子 能 科学 技术 ,2015(9):1624-1628. DOT:10.7538/yzk.2015.49.09.1624.) 


[6] Zeguang Li, Jun Sun, MaLin, et al.Nuclear engineering and design,2020,361(May):110569.1-110569.12. 
DOI:10.1016/j.nucengdes.2020.110569. 

[7] Jian Li, Qiu Zhou, Jian Mou, et al.Annals of nuclear energy,2020,142(Jul.): 107382. 1-107382.8. DOI:10.1016/j.anucene.2020.107382. 
[8] LI Ting, ZHUANG Kun, SHANG Wen, et al. Nuclear Techniques,2020,43(8):34-42(in Chinese). 
DOTI:10.11889/j.0253-3219.2020.hjs.43.080006. 

( 李 婷 , 庄 坤 ,尚文 ,等 . 核 技术 ,2020,43(8):34-42. DOI:10.11889/j.0253-3219.2020.hjs.43.080006.) 


[9] JIANG Wei, ZHANG Lu, YU Rui, et al.Atomic Energy Science and Technology,2022,56(5):961-968(in Chinese). 
DOI:10.7538/yzk.2022.youxian.0173. 
( 姜 韦 , 张 璐 ,于 锐 ,等 .原子 能 科学 技术 ,2022,56(5):961-968. DOI:10.7538/yzk.2022. youxian.0173.) 


[10] Anatoli Tsiboulia, Oleg A. Elovsky.NEA/NSC/DOC/(95)03/IL 1987. 


[11] Chadwick M B, Herman M, Oblo Insky P, et al. Nuclear Data Sheets, 2011, 112(12): 2887—2996. 


— 


2] Zhigang Ge, Ruirui Xu, Haicheng Wu, et al. EPJ Web of Conferences, 2020, 239 : 09001-. 
[13] Brown D A, Chadwick M B, Capote R, et al. Nuclear Data Sheets, 2018, 148: 1-142. 


[14] A.J.M. Plompen, O. Cabellos, C. De Saint Jean, et al. Eur. Phys. J. A. 56 (2020) 181, 


https://doi.org/10.1140/epja/s10050-020-00141-9. 


Study on criticality benchmark testing in spatial reactors 


based on different evaluated nuclear databases 


Chen Xiang! Ma Xubo! Fu Liantao! Ma Xudong! 
(1.North China Electric Power University, School of Nuclear Science and Engineering, Beijing 102206) 


Abstract: In order to improve the reliability of neutronics simulation for space reactors, based on the criticality 
benchmark experiment for space reactors in the International Criticality Safety Handbook (HMF005 benchmark 
experiment), the Monte Carlo program is used to model the core arrangement in the HMF005 benchmark 
experiment, and the deviation of the ke calculated values from the experimental values under six core 
arrangements is calculated by combining the different evaluated nuclear databases and the ratios between the 
calculated and experimental values are analyzed in depth. mean, standard deviation, y? and root mean square error 
are analyzed in depth. It is found that the deviation of the core ke results from the experimental values calculated 
by the four types of evaluation nuclear databases, ENDF/B-VII.1, CENDL-3.2, ENDF/B-VIII.0 and JEFF-3.3, is 
within the triple uncertainty, with the calculated results of ENDF/B-VII.1 being the closest to the experimental 
values, and the deviation of CENDL-3.2 from the experiments being a little bit larger. The deviation of the 
calculated results from the experimental values of the CENDL-3.2 library is also shown in the table below. 
CENDL-3.2 library calculations deviate slightly from the experimental values mainly due to the molybdenum and 
beryllium nuclides in the evaluation library, and the isotope of molybdenum has the largest effect of ??Mo. 
Sensitivity analysis was also carried out, and the total reaction sensitivity coefficients of ?Mo and ?Be were mainly 
related to the elastic scattering sensitivity coefficient, and the difference in radiation capture sensitivity between 
CENDL-3.2 and ?Be in ENDF/B-VII.1 was large, and the deviation between the two could be up to 34.9%. 

Key words: critical calculations; effective multiplication factor; evaluation of the nuclear database; Monte Carlo 


simulation 


! Received: 2023-12-29 Modified: 2024-3-7 

Foundation item: The National Natural Science Foundation of China (11875128) 

First author: Chen Xiang (1999-), Male, Tujia Family, From Shizhu, Chongqing, Master's Degree. Research Interests: Nuclear Energy 
Science and Engineering. E-mial:120222212004@ncepu.edu.cn 

Corresponding author: Ma Xubo maxb@ncepu.edu.cn 


